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Abstract: We analyse the effective potential of the scalar wave equation near generic 
space-time singularities of power-law type (Szekeres-Iyer metrics) and show that the ef- 
fective potential exhibits a universal and scale invariant leading inverse square behaviour 
~ in the "tortoise coordinate" x provided that the metrics satisfy the strict Dominant 
Energy Condition (DEC). This result parallels that obtained in Q for probes consisting 
of families of massless particles (null geodesic deviation, a.k.a. the Penrose Limit). The 
detailed properties of the scalar wave operator depend sensitively on the numerical coeffi- 
cient of the x~^-term, and as one application we show that timelike singularities satisfying 
the DEC are quantum mechanically singular in the sense of the Horowitz-Marolf (essential 
self-adjointness) criterion. We also comment on some related issues like the near-singularity 
behaviour of the scalar fields permitted by the Priedrichs extension. 
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1. Introduction 

The study of scalar field propagation in non-trivial curved (and possibly singular) back- 
grounds is of fundamental importance in a variety of contexts including quantum field 
theory in curved backgrounds, cosmology, the stability and quasi-normal mode analysis of 
black hole metrics etc. 

Typically, this is studied within the context of a particular metric or class of metrics. 
For certain purposes, however, only the knowledge of the leading behaviour of the metric 
near a horizon or the singularity is required. In that case, one can attempt to work with 
a general parametrisation of the metric near that locus and, in this way, ascertain which 
features of the results that have been obtained previously for particular metrics are special 
features of those metrics or valid more generally. 

In particular, practically all explicitly known metrics with singularities are of "power- 
law type" in a neighbourhood of the singularity (instead of showing, say, some non- 
analytic behaviour). In the spherically symmetric case, the leading behaviour of generic 
metrics with such singularities of power-law type is captured by the 2-parameter family 

ds"^ = r]xP{-dx^ + dy^) + x'^dQ^ (1.1) 
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of Szekeres-Iyer metrics ||2|, ^, ^. The singularity, located in these coordinates at x = 0, is 
timelike for rj = —1 and spacelike for rj = +1. This class of metrics thus provides an ideal 
laboratory for investigating the behaviour of particles, fields, strings, ... in the vicinity of 
a generic singularity of this type. 

A first investigation along these lines was performed in |l| in the context of the 
Penrose Limit, i.e. of probing a space-time via the geodesic deviation of families of massless 
particles. There it was shown that the plane wave Penrose limits, 

ds"^ = g^ydx^dx" 2dudv + Aab{u)x''x^du^ + dx'^ , (1.2) 

of metrics with singularities of power-law type have a universal ii~^-behaviour near the 
singularity, Aab{u) ~ u~'^, provided that the near-singularity stress-energy (Einstein) tensor 
satisfies the strict dominant energy condition (DEC). This behaviour, which is precisely 
such that it renders the plane wave metric scale invariant Q , had previously been observed 
in various particular examples and is thus now understood to be a general feature of this 
class of singularities. 

It is then natural to wonder whether a similar universality result can be established in 
other circumstances or for other kinds of probes and if, analogously, some energy condition 
plays a role in establishing this. If one considers e.g. the Klein-Gordon equation = for 
scalar fields, it is not difficult to see 0, |8| that under certain conditions the scalar effective 
potential V^s for general metrics with singularities of power-law type displays an inverse 
square behaviour, Ves{x) ~ x~'^, near the singularity. This observation was then used in 
Q to study the quasi-normal modes for black holes with generic singularities of this type. 

The purpose of this note is to study other aspects and consequences of this universality. 
In particular, we will first show that the results obtained in [l[ , namely the scale invariant 
inverse square behaviour of the wave profile Aab{u), as well as a crucial |9|, |^ lower bound 
on the coefficients, have a precise and rather striking analogue in the case of a scalar field. 
Schematically, this analogy can be expressed as 

^ ■ ^ ^ J ^ab{u) CaSabU'"^ scale invariant (ca > -1/4) 

strict DEC \ t \ -9 1 / (1-3) 

I Vcq{x) — > cx scale mvariant (c > —1/4) 

Once again this shows that this inverse square behaviour, that had been observed before 
in various specific examples in a variety of contexts, is a general feature of a large class 
of space-time singularities. The precise statements are derived in sections 2.2 and 2.3 and 
discussed in section 2.4, while sections 2.5 and 2.6 deal with minor variations of this theme. 

We hasten to add that if such an inverse square behaviour were universally true without 
any further qualifications then it would probably have to be true on rather trivial (dimen- 
sional) grounds alone. What makes the results obtained here and in Q more interesting 
is that a priori in either case a more singular behaviour can and does occur and is only 
excluded provided that some further (e.g. positive energy) condition is imposed. 

The significance of the x~^-behaviour is that (as anticipated in (|1.3D ), the correspond- 
ing Schrodinger operator —d"^ + cx~'^, to which we will have reduced the Klein-Gordon 
operator, defines a scale invariant (c is dimensionless) "conformal quantum mechanics" 
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flOfl problem. Thus, here and in we find a rather surprising emergence of scale invari- 
ance in the near-singularity limit. One implication of this scale invariance in the plane wave 
case, discussed in [11 1, is that it leads to a Hagedorn-like behaviour of string theory in this 
class of backgrounds that is quite distinct from that in plane wave backgrounds with, say, 
a constant profile and more akin to that in Minkowski space. It would be interesting to 
explore other consequences of this near-singularity scale invariance. 

This class of scale invariant models has recently also appeared and been discussed in 
various other related settings, most notably in the analysis of the near- horizon (rather 
than the near-singularity) properties of black holes, see e.g. |12, 14, 15, 1£], where the 
emergence of scale invariance can largely be attributed to the near-horizon AdS geometry, 
as well as in quantum cosmology |17]. 

Having reduced the Klein-Gordon operator to the Schrodinger operator —9^ -|- crc~^ 
(after a separation of variables and a unitary transformation), one can then turn to a 
more detailed spectroscopy of the Szekeres-Iyer metrics by analysing the properties of this 
operator. Indeed, as is well known, the inverse square potential is a critical borderline case 
in the sense that the functional analytic properties of this operator depend in a delicate way 
on the numerical value of the coefficient c. This value, in turn, depends on the dimension 
d (number of transverse dimensions) and the Szekeres-Iyer parameter q (it turns out to be 
independent of p, while the corresponding coefficients Ca in the Penrose limit case typically 
depend on {p,q) and d). 

As one application, we will analyse the Horowitz-Marolf criterion |18| for general singu- 
larities of power-law type. Horowitz and Marolf defined a static space-time to be quantum 
mechanically non-singular (with respect to a certain class of test fields) if the evolution of 
a probe wave packet is uniquely determined by the initial wave packet (as would be the 
case in a globally hyperbolic space-time) without having to specify boundary conditions 
at the classical singularity. This criterion can be rephrased as the condition that the (spa- 
tial part of the) Klein-Gordon operator be essentially self-adjoint (and thus have a unique 
self-adjoint extension). 

While such a necessarily only semi-classical analysis is certainly not a substitute for a 
full quantum gravitational analysis, it nevertheless has its virtues since one can learn what 
kind of problems persist, can arise or can be resolved when passing from test particles to 
test fields. 

Intuitively one might expect a classical singularity with a sufficiently "positive" (in an 
appropriate sense) matter content to remain singular even when probed by non-stringy test 
objects other than classical point particles. This line of thought was one of the motivations 
for analysing the Horowitz-Marolf criterion in this framework, and we will indeed be able 
to show (section 3.4) that 



metrics with timelike singularities of power-law type satisfying the strict Dom- 
inant Energy Condition remain singular when probed with scalar waves. 

A second issue we will briefly address is that of the allowed near-singularity behaviour of 
the scalar fields for a given self-adjoint extension (section 3.5). A priori, one might perhaps 
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expect a sufficiently repulsive singularity to be regular in the Horowitz-Marolf sense simply 
because the corresponding unique self-adjoint extension forces the scalar field to be zero 
at the singularity, thus in a sense again excluding the singularity from the space-time. It 
is also possible, however, and potentially more interesting, to have a self-adjoint extension 
with scalar fields that actually probe the singularity in the sense that they are allowed to 
take on non-zero values there. We propose to call such singularities "hospitable" , establish 
once again a relation, albeit not a strict correlation, with the DEC, and show among other 
things that, in a suitable sense, half of the Horowitz-Marolf regular power-law singularities 
are hospitable whereas the others are not. 



2. Universality of the Effective Scalar Potential for Power-Law Singulari- 
ties 

2.1 Geometric Set-Up 

Even though we will ultimately be interested in the properties of the scalar wave (Klein- 
Gordon) equation (□ — m^)(/) = in the Szekeres-Iyer metrics (|l . l|) , to set the stage it 
will be convenient to begin the discussion in the more general setting of metrics with 
a hypersurface orthogonal Killing vector. The general set-up here and in section 3.1 is 
modelled on the approach of [|l9| (with minor adaptations to allow for both timelike and 
spacelike singularities). 

We begin with the n-dimensional metric 

ds^ = rja^dy^ + hijdx^dx^ (2.1) 

where a and hij are independent of y, = dy is a hypersurface orthogonal Killing vector 
with norm ^^^^ = rja^, and thus timelike (spacelike) for t] = —1 {t] = +1). Correspondingly 
we assume that the metric hij induced on the hypersurfaces = S of constant y is 
Riemannian (Lorentzian) for r] = —1 [r] = +1). 

Denoting the covariant derivatives with respect to the metric hij by Di, the wave 
operator 

□ ^ J^ d^./^d^gf^'^d, (2.2) 
a/— det g 

is easily seen to take the form 

□ = a-^(ry5^ + aD'aDi) . (2.3) 
Thus the massive wave equation (□ — m^)0 = can be written as 

d^y4> = , (2.4) 

where A is the operator 

A = rjaD'aDi - T]a^Tr? . (2.5) 
Assuming now spherical symmetry, the metric takes the warped product form 

ds^ = rja{xfdy^ - r]b{xfdx^ + cixfdnl (2.6) 
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where dO^, d = n — 2, denotes the standard metric on the d-sphere S"^. It wih be apparent 
from the fohowing that the assumption of spherical symmetry could be relaxed - we will 
only use the warped product form of the metric in an essential way. 

We could fix the residual x-reparametrisation invariance by introducing the "area ra- 
dius" r = c(x) as a new coordinate. However, for the following it will be more convenient 
to choose the gauge a{x) = b{x) (i.e. x is a "tortoise coordinate" for rj = —1 respectively 
"conformal time" for rj = +1), 

ds"^ = j]a{xf{-dx^ + dy^) + c{xfdill . (2.7) 

Then the operator A is 

A = -a'^dxadx + 'na^c~'^Ad - rja^m^ , (2.8) 

where a{x) = c(x)°' and denotes the Laplacian on S'^. 

To put A into standard Schrodinger form, we transform from the functions (j){x) to 
the half-densities (cf. ( p.Sp ) 0(x) = a^^'^(j){x). The corresponding unitarily transformed 
operator A is 

A = a^/^Aa-^/^ = -dl + V + va^c'^Aa - ria^m^ 
V{x) = a{xr^/\dla{xf/^) . (2.9) 

After the usual separation of variables in the y-direction, 

<^(y,x,r) = e-^^y^(x,r) , (2.10) 

and the decomposition into angular spherical harmonics Y£^{9"'), with 

-Arfy,^(r) = ilY.rnie") 

£l = £{£ + d-l) , (2.11) 

the Klein-Gordon equation for the metric ( p.7| ) reduces to a standard one-dimensional 
time-independent Schrodinger equation 

[-dl + Kff,£(x)] <^(x) = E^{x) (2.12) 

{(j){x) = (l)E,e,in{x)) with effective scalar potential 

Ves,i{x) = V{x)-r]a{xf{£jc{x)-^ + m^) . (2.13) 

2.2 The Effective Scalar Potential for Power-Law Singularities 

The leading behaviour of generic (spherically symmetric) metrics with singularities of 
power-law type^, i.e. metrics of the general form 

ds^ = -dt^ + [t- r(r)]2V(r, ifdr"^ + [t- T{r)f^g{r, tfd^l , (2.14) 

^Such metrics encompass practically all explicitly known singular spherically symmetric solutions of 
the Einstein equations like the Lemaitre-Tolman-Bondi dust solutions, cosmological singularities of the 
Lifshitz-Khalatnikov type, etc. On the other hand, this class of metrics does prominently not include the 
BKL metrics describing the chaotic oscillatory approach to a spacelike singularity. Whether or not such 
a behaviour occurs depends in a delicate way on the matter content, see e.g. |2l| and references therein. 
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with / and g functions of r and t that are regular and non-vanishing at the location t = r(r) 
of the singularity, is captured by the 2-parameter family of Szekeres-Iyer metrics |2|, ^] (see 
also HI and the generalisation to string theory backgrounds discussed in Q) 

ds'^ = r]xP{-dx'^ + dy^) + x'^dQ^ . (2.15) 

The Kasner-like exponents p,q G M. characterise the behaviour of the geometry near the 
singularity at 2; = 0. This singularity is timelike for t] = —1 (x is a radial coordinate) 
and spacelike for ry = -|-1 (with x a time coordinate). In particular, these metrics possess 
the hypersurface orthogonal Killing vector dy, and are already in the "tortoise" form (p.7|), 
with a(x)^ = x^ and c(x)^ = x'^. Thus we can directly read off the effective scalar potential 
from the results of the previous section. 

From (p.g|), we deduce, with o"(x) = x'^''/^, that 

V{x) = s{s-l)x-^ « = X ■ ^^-^^^ 

Thus, from ( |2.13| ) we find (see also ^) 

Ves,i{x) = s{s - l)x-^ - v^lxP-'^ - rjm^xP (2.17) 

We are interested in the leading behaviour of this potential as 2; ^ (subdominant terms 
can in any case not be trusted as we have only kept the leading terms in the metric ( |2.15| )). 
For the time being we will consider the massless case = (see section 2.5 for 7^ 0). 

Provided that s(s — 1) 7^ 0, which term in ( |2.17| ) dominates depends on p and q. When 
q < p + 2, one finds 

q<p + 2: VesAx)^s{s-l)x-^ . (2.18) 

The two salient features of this potential are the inverse square behaviour and a coefficient 
c that is bounded from below by —1/4, 

c = s{s-l)>-^ , (2.19) 

with equality for s = 1/2, i.e. q = 2/d. 

As mentioned in the introduction, the significance of the x^^-behaviour is that it defines 
a scale invariant "conformal quantum mechanics" jl^ problem, discussed more recently in 
related contexts e.g. in |12, 14, |l^, 16, |l^. Moreover, for practical purposes |8|, |2^ the 
virtue of the (as opposed to a more singular) behaviour is that it leads to a standard 
regular-singular differential operator. 

The significance of the bound on c is that in this range the operator —5^ -|- c/x^ is 
positive, as can be seen by writing 

-dl + s{s - l)x-2 = {d^ + sx~^){-d^ + sx"^) = {-d^^ + sx-^)^{-d^ + sx"^) . (2.20) 

When q = p + 2, the metric is conformally fiat, both terms in ( |2.17| ) contribute equally, 
and one again finds the x~^-behaviour 

q = p + 2: V^sM^cx-^ , (2.21) 
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where now 

c = s{s-l)-r^tl . (2.22) 

Thus in this case c is still bounded by —1/4 for timelike singularities, while c can become 
arbitrarily negative for sufficiently large values of for r] = +1. 

Once q > p + 2, the second term in ( 2.171 ) dominates (for £^ / 0), and one finds the 



more singular leading behaviour 

q>p + 2: V^Q^iix) -rji^ x'^"" a > . (2.23) 

Examples of metrics with q < p + 2 are the Schwarzschild and Friedmann- Robertson- 
Walker (FRW) metrics and indeed, as we will recall below, all metrics satisfying the strict 
Dominant Energy Condition. 

In particular, for the (d + 2)-dimensional (positive or negative mass) Schwarzschild 
metric, one has 

Schwarzschild : p = — - — Q = , (2.24) 

as is readily seen by expanding the metric near the singularity and going to tortoise co- 
ordinates. Thus the Schwarzschild metric has s = 1/2 and c takes on the d-independent 
extremal value c = —1/4, as observed before e.g. in ||2^ , ^ in related contexts. 

For decelerating cosmological FRW metrics, with cosmological scale factor (in comov- 
mg time) t'', < h < I, 

h = — ^ , (2.25) 

{d + l){l + w) ' ^ ^ 

with w the equation of state parameter, P = wp, one finds ||5|, |l| 

2h 

FRW: p = q = -^ , (2.26) 

as can be seen by going to conformal time. A routine calculation shows that the above 
result (|2.18 ) for the purely x-dependent part of the effective potential (with x conformal 



time) is actually an exact result, and not an artefact of the near-singularity Szekeres-Iyer 
approximation. 

It remains to discuss the case when q < p + 2, so that the first term in ( |2.17| ) would 
be dominant, but the coefficient s{s — 1) = 0. When s = 0, then one has q = and this 
is generally interpreted |Q] as corresponding not to a true central singularity (as the radius 
of the transverse sphere remains constant as x — > 0) but as a shell crossing singularity. 

The other possibility is s = 1, i.e. q = A/d. This is a case in which (because of the 
cancellation of the leading terms) subleading corrections to the metric ( p. 15 ) can become 



relevant and should be retained. An example of metrics with s = 1 is provided by FRW 
metrics with a radiative equation of state. Using ( 2.26| ), one has 



which is precisely the equation of state parameter for radiation. However, as follows from 
the remark above, in this special case the vanishing of the effective potential for p = q \s 
actually an exact result. 
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2.3 The Significance of the (Strict) Dominant Energy Condition 

We have seen that genericaUy the leading behaviour of the scalar effective potential near 
a singularity of power-law type is either ~ x~'^ or ~ x:P~'^. We will now recall from |l| 
that the latter behaviour can arise only for metrics violating the strict Dominant Energy 
Condition (DEC). While there is nothing particularly sacrosanct about the DEC, and 
other energy conditions could be considered, the DEC appears to play a privileged role in 
exploring and understanding the (p, q)-plane of Szekeres-Iyer metrics. 

The Dominant Energy Condition on the stress-energy tensor (or Einstein tensor 
G^u) [^] requires that for every timelike vector v^, T^j^yV^v'^ > 0, and T^v'^ be a non- 
spacelike vector. This may be interpreted as saying that for any observer the local energy 
density is non-negative and the energy flux causal. 

The Einstein tensor of Szekeres-Iyer metrics is diagonal, hence so is the corresponding 
stress-energy tensor. In this case, the DEC reduces to 

P > \Pi\ , (2.28) 

where —p and Pi, i = 1, . . . ,d + 1 are the timelike and spacelike eigenvalues of Tu re- 
spectively. We say that the strict DEC is satisfied if these are strict inequalities and we 
will say that the matter content (or equation of state) is "extremal" if at least one of the 
inequalities is saturated. 

Now it follows from the explicit expression for the components 

G^^ = -\d{d - l)x-'i - \r]dq{{d - l)q + 2p)x-(P+2) 

Gl = -\d{d - l)x-« + l^ndqilp + 4 - (d + l)g)x-(P+2) (2.29) 

of the Einstein tensor that for q> p + 2 the relation between —p and the radial pressure Pr 
(identified with G% and - which is which depends on the sign of rj) becomes extremal 
asx^O |,|, 

q>p + 2: Gl-Gl^Q ^ p + Pr^O . (2.30) 

Put differently, q < p + 2 is a necessary condition for the strict DEC to hold, and thus for 
metrics satisfying the strict DEC the leading behaviour of the effective potential is always 

As an aside, we note that it follows from ( |2.29| ) that precisely those metrics that satisfy 
the physically more reasonable (non-negative pressure) and more common extremal near- 
singularity equation of state p = +Pr have q = 2/d, i.e. s = 1/2, leading to the critical 
value c = —1/4 frequently found in applications (to e.g. Schwarzschild-like geometries). 

2.4 Comparison with Massless Point Particle Probes (the Penrose Limit) 

In the previous section we have established that 

1. for metrics with singularities of power-law type satisfying the strict DEC the leading 
behaviour of the scalar effective potential near the singularity is 

Veff,^(x) ^ (2.31) 
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2. this class of potentials is singled out by its scale invariance; 

3. the corresponding coefficient c of the effective potential is bounded from below by 
— 1/4 unless one is on the border to an extremal equation of state. 

These observations bear a striking resemblance to the results obtained recently in ||T| 
in the study of plane wave Penrose limits 

ds'^ = g^^dxf'dx" 2dudv + Aab{u)x''x''du'^ + dx^ , (2.32) 

of space-time singularities. Namely, it was shown in [P that 

1. Penrose limits of metrics with singularities of power-law type show a universal Mis- 
behaviour near the singularity, 

Aab{u) CaSabU'"^ , (2.33) 



provided that the strict DEC is satisfied; 



such plane waves are singled out Q by their scale invariance, reflected e.g. in the 
isometry {u,v) {\u,X~^v) of the metric ( |2.32| , 2.33| ); 

the coefficients Ca (related to the harmonic oscillator frequency squares by Ca = — w^) 
are bounded from below by —1/4 unless one is on the border to an extremal equation 
of state. s 



The similarity of these two sets of statements is quite remarkable because the objects 
these statements are made about are rather different. For example, the potential is that of 
a one-dimensional motion on the half line in one case, and that of a d-dimensional harmonic 
oscillator (with time-dependent frequencies) in the other. 

The analogy with the above statements about scalar effective potentials is brought 
out even more clearly if one reinterprets |5|, |^ the Penrose limit in terms of null geodesic 
deviation in the original space-time. Then this result can be rephrased as the statement that 
the leading behaviour of the geometry as probed by a family of massless point particles 
near a singularity is that of a plane wave with a u~'^ geodesic effective potential. The 
analogy with the results of the previous section should now be apparent. 

One minor difference between the results obtained here and those of 0] is that in the 
case of Penrose limits the strict DEC needed to be invoked only in the case of spacelike 
singularities, rj = +1, timelike singularities always giving rise to plane waves with a Mis- 
behaviour. This should be regarded as an indication (cf. the discussion in [||, Section 4.4]) 
that scalar waves are better probes of timelike singularities than massless point particles. 

^One significance of this bound on the Ca is that in this range one can consider the possibility to extend 
the string modes across the singularity at m = M. 
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2.5 Massive Scalar Fields and Geodesic Incompleteness 

The simple above analysis can evidently be generalised in various ways, e.g. by considering 
other kinds of probes. We will briefly comment on the two most immediate generalisations, 
namely massive and non-minimally coupled scalar fields. 

We begin with a massive scalar for which the effective potential is 

Ves/ix) = s{s - - r/^rfX^-'? - rjm^xP (2.34) 

For the mass term to be relevant (dominant) as x ^ it is clearly necessary that p < —2 
and q < 0. Intuitively one might expect a mass term to be irrelevant at short distances 
near a singularity. This expectation is indeed borne out: as we will now show, for metrics 
satisfying the above inequalities the would-be singularity at x = is at infinite affine 
distance for causal geodesies so that such space-times are actually causally geodesically 
complete. 

Null geodesies were analysed in [p. Here we generalise this to causal geodesies. In 
terms of the conserved angular and y-momentum L and P, the geodesic equation for the 
metric ( p. 15 ) reduces to 



X 



^ - P^x-^P + r]L'^x-P~'^ + r]ex-P , (2.35) 



where e = (e = 1) for null (timelike) geodesies respectively. 

For T] = —1, if the first term in (|2.35| ) is sub-dominant the geodesic effective potential is 
repulsive (e.g. via the angular momentum barrier) and the geodesies will not reach x = 0. 
Thus generic timelike geodesies will reach x = only if (p, q) lie in the positive wedge 
bounded by the lines p = and p = q. Radial null geodesies do not feel any repulsive force, 
and solving 

^ x-^P ^ x{u) ^\ ^ ^ (2.36) 

I exp u p = —1 

shows that x = is reached at a finite value of the affine parameter only for p > —1. We 
thus conclude that Szekeres-Iyer metrics with ry = — 1 and p < —1 are causally geodesically 
complete. In particular, therefore, the mass term in the scalar effective potential is sub- 
dominant for metrics with honest timelike power-law singularities, and for all such metrics 
the scalar effective potential has the same leading behaviour as in the massless case. 

For 7] = +1, the situation is more complex as all three terms in ( 2.35| ) are positive. If 



the first term dominates, either because of suitable inequalities satisfied by {p, q) or, for 
any (p, q), because one is considering radial null geodesies, the analysis and conclusions are 
identical to the above. Namely, x = is at finite affine distance for p > —1. Analogously, if 
the second term dominates (e.g. for angular null geodesies) one finds the condition p + q > 
—2, and if the third term dominates one has p > —2. Since one needs p < —2 for the mass 
term to dominate in the scalar effective potential, only the second case is possible. But 
then the condition p + q > —2, with p < —2, implies q > 0, so that the angular momentum 
term in the effective potential dominates the mass term. 

We thus conclude that, for both r] = +1 and ij = —1, the mass term is always sub- 
dominant for metrics that are causally geodesically incomplete at x = 0. 
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As an aside we note that the Szekeres-Iyer metrics for which the mass term does 
dominate {p < —2 and g < 0), in addition to being non-singular, also necessarily violate 
the strict DEC. 



2.6 Non-Minimally Coupled Scalar Fields 

We will now very briefly also consider a non-minimally coupled scalar field 

(□ - = . (2.37) 

The Ricci scalar of the Szekeres-Iyer metric ( 2.15| ) is 

R = d{d - l)x-'^ - lr]{Ap + 4qd - d{d + l)q'^)x-'-P+^^ , (2.38) 

where once again only the leading order term should be trusted and retained. Thus the 
new effective potential 

V^s/x) = VesA^) - r]CxPR (2.39) 

is again a sum of two terms, proportional to and x^""^ respectively, so that the dominant 
behaviour is still ~ x~'^ provided that the metric does not violate the strict DEC. For 
q < p + 2 and the conformally invariant coupling 

one finds 

Note that with this conformally invariant coupling the coefficient c now depends on p — 
q rather than on q. The appearance of (p — q) could have been anticipated since for 
p = q the Szekeres-Iyer metric is conformal to an x-independent metric, and hence a 
conformal coupling cannot generate an x-dependent effective potential. Note also that 
for the conformal coupling (and, indeed, generic values of ^) the coefficient c is no longer 
bounded by —1/4 so that the Schrodinger operator is no longer necessarily bounded from 
below. 



3. Self- Adjoint Physics of Power-Law Singularities 

In the previous section we have determined the leading behaviour of the scalar wave op- 
erators near a power-law singularity. In this section we will now study various aspects of 
these operators. 

3.1 Functional Analysis Set-Up 

In order to analyse the properties of the wave operator, we will need to equip the space of 
scalar fields with a Hilbert space structure. We will be pragmatic about this and introduce 
the minimum amount of structure necessary to be able to say anything of substance. 
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We thus return to the discussion of section 2.1, now being more specific about the 



spaces of functions the various operators appearing there act on [19|, beginning with the 
operator A introduced in (^), 

A = r]aD'-aDi-r]a'^m^ . (3.1) 

Since D^Di is symmetric (formally self-adjoint) with respect to the natural spatial density 
\/—r] det h induced on the slices S of constant y by the metric ( |2.1| ), the operator A is 
symmetric with respect to the scalar product 

(</>!, <^2) = I d"-'xacl)lct>2 

o- = a"^y^^Miet7i = T/y^det^g^J' , (3.2) 

on D{A) = Co-(S), 

(yl0i,02) = (0i,A02) . (3.3) 
Moreover, for 7] = —1 the operator A is positive, 

r] = -l =^ (0, A0) > . (3.4) 

We are thus led to introduce the Hilbert space L^(S, crd^^^x) of functions on S square 
integrable with respect to the above scalar product. 

Passing to spherically symmetric metrics ( |2.6| ) in the tortoise gauge (^), A takes the 
form (|2]|) 

A = —a~^dxcrdx + rja^c^'^Aci — rja^m^ , (3.5) 



where a{x) = c{x) . Since A is symmetric with respect to the scalar product (p. 21) , the 
unitarily transformed operator 

A = a'/^Aa-'/' , (3.6) 

acting on the half-densities 

0(x) = a{x)^/^(P{x) , (3.7) 
is symmetric with respect to the corresponding "flat" {cr{x) 1) scalar product 

< 01,(^2 >:= / dxdn (f>l4>2 = {(1)1,(1)2) ■ (3.8) 



We now assume that the metric develops a singularity at some value of x, where e.g. the 
area radius goes to zero, r = c{x) — >■ 0, which we may as well choose to happen at a; = 0. 
Thus we consider x € (0,cxd) and take S = R"'~^\{0}, parametrised by x and the angular 
coordinates. 

Then the initial domain of i is D{A) = C^(M"-i\{0}) or D{A) = C^(M+) ® C°°(S"^), 
which are dense in the unitarily transformed Hilbert space 

l2(M"-1\{0}, dx dn) ^ L\R+,dx) 8) L\S'^, dQ) . (3.9) 
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Decomposing the second factor into eigenspaces of the Laplacian A^^ on S"^, 

oo 

L'^{R+,dx)0L'^{S'^,dn) = ^Li , (3.10) 
and defining = D D Li, one has 

A\^^=Ae(g)I , (3.11) 

where 

Ae = -dl + Kff/(x) (3.12) 

with 14fr/(x) given in ( |2.13| ). 

Questions about the original operator A can thus be reduced to questions about the 
family {^f} of standard Schrodinger-type operators. For example, to show that A is 
essentially self-adjoint on D{^A) it is sufficient to prove that, for each ^, A^ is essentially 
self-adjoint on C^(M+). 

While one can analyse this question of self-adjointness just as readily for r] = +1 as for 
r] = —1, the physical significance of this condition in the case of spacelike singularities is 
not clear to us. Thus we will focus on static space-times with timelike singularities in the 
following and set rj = —1. An extension of the general formalism to stationary non-static 
space-times is developed in |24|. 



We conclude this section with a comment on the choice of Hilbert space structure. 
The Hilbert space introduced above is certainly a natural choice, but not the only one 
possible. Based on physical requirements such as the finiteness of the energy of scalar field 
probes, other (Sobolev) Hilbert space structures have been proposed in the literature - see 
e.g. |2^, 26 1 . The energy is, by definition, 

E[<))] = [ Vkd^'-^x T^M)en'' , (3.13) 

where T^^{(p) is the stress energy tensor of the scalar field, = dy is the timelike Killing 
vector, and n the unit normal to S. In the present case this reduces to 

E[^] = j ad^-^x Tyy , (3.14) 

which identifies Tyy as the energy density with respect to the measure ad^^^x employed 
above |2(:]. For a minimally coupled complex scalar field one has 

Tyy = ^[^yCj)*^yCj) + a^h'^^^C^*^jC|)\ . (3.15) 

Thus, with an integration by parts (certainly allowed for (p € D{A)) the energy can be 
written as 

E[^] = [ ad'^-^x {dy(i)*dy^ + <t)*A4)) 

= {dy(t),dy(t)) + {(t),A(l)) . (3.16) 
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For a comparison of the two definitions ( 3.14| ) and ( ^.16 ) of the energy and the role of bound- 



ary terms, see e.g. the discussion in [27| and the comment in section 3.5 below. Adopting 



the expression ( 3.16| ) as the definition of the energy suggests introducing a Sobolev struc- 



ture on the space of scalar fields using the quadratic form 

Qa{^) = {(P,A^) (3.17) 



associated to A, via |23, 26 1 



U\\l. = {(l>A) + QA{<t>) , (3.18) 

thus enforcing the condition that the energy be finite. For present purposes we simply note 
that at least for the Friedrichs extension Ap of A, based on the closure of the quadratic 
form QAiff) with respect to the norm, the resulting potential energy QAp{(k) is finite 
(and positive) by definition without having to invoke Sobolev spaces (see also the discussion 
in p8| , pS|]).'^ We will use specifically this extension in the discussion of section 3.5 below. 

3.2 Essential Self-Adjointness and the Horowitz-Marolf Criterion 

The spatial part A of the wave operator is real and symmetric (with respect to an ap- 
propriate scalar product on a domain of A)^ and as such has self-adjoint extensions, 
each leading to a well defined (and reasonable |^8|) time-evolution. If the self-adjoint ex- 
tension is not unique, however, i.e. if the operator is not essentially self-adjoint, then also 
the corresponding time-evolution is not uniquely determined. Thus the Horowitz-Marolf 
criterion (unique time-evolution without having to impose boundary conditions at the 
singularity) amounts to the condition that the operator A be essentially self-adjoint. 

To test for essential self-adjointness pH], one can e.g. use [18| the standard method of 



Neumann deficiency indices or the Weyl limit point - limit circle criterion (employed in 
this context in [0]). Roughly speaking, in order for A to be essentially self-adjoint the 
(effective) potential V'eff,^ appearing in the operator A^ has to be sufficiently repulsive near 
X = to prevent the waves cj) from leaking into the singularity. 

Concretely, in the present case, where we only have control over the operator A near 
the singularity at x = 0, the criteria for the operator to be essentially self adjoint 
on C^(M+) at x = boil down to the following elementary conditions on the effective 
potential = Kg ^ IsO] : 



If 



W{x) > ^x-^ (3.19) 



near zero, then —5^ -|- W{x) is essentially self-adjoint at x = 0. 
If for some e > 



W{x) < (^-e]x~^ (3.20) 



■^Working with such a Sobolev space structure is certainly possible but also complicates the determination 
of self-adjoint extensions of A, since e.g. studying the closure of A now involves studying the sixth order 
operator A^, arising from the term — {Acl>,A(f>) + {A(j),A?(j)) in the operator norm. 
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(in particular also if is decreasing) near x = 0, then —d'^+W{x) is not essentially 
self- adjoint at x = 0. 



The significance of the factor 3/4 can be appreciated by looking at the critical (and 
relevant for us) case of an inverse square potential 

W{x) = s{s-l)x-^ . (3.21) 
In this case the leading behaviour of the two linearly independent solutions of the equation 

{-dl + W{x))Mx) = ^Mx) (3.22) 
near x = is given by the two linearly independent solutions of the equation 

{-dl + W{x))M^)=0 , (3.23) 

namely 

4)0 ~ x' or 4>o ~ x^-' (3.24) 
Thus both solutions are square integrable near x = when 2s > —1 and 2(1 — s) > —1, or 

<s< ^ ^ s{s-l)<^ , (3.25) 

and in this range of c = s(s — 1) the potential is limit circle and the self-adjoint extension 
is not unique. Conversely, it follows that for c > 3/4 the solutions of equation (|3.22| ) for 
X = ±1 (which are necessarily complex linear combinations of the two linearly independent 
real solutions) are not square-integrable near x = 0. Thus the deficiency indices are zero 
and the operator is essentially self-adjoint for c > 3/4. 

Even when there are two normalisable solutions, all is not lost however, as it may 
be indicative of the possibility (or even necessity) to continue the fields and/or the metric 



through the singularity [ 22 1 . Evidently, such an analytic continuation requires some thought 
(to say the least) in the case of Szekeres-Iyer metrics with generic (non-rational) values of 
p and q. 

3.3 The Horowitz-Marolf Criterion for Power-Law Singularities 

In the case at hand, timelike singularities of power-law type, the effective potential is given 



by (|2.17D with rj = —1 and s = qd/A. We had already seen in section 2.5 that the mass 
term is never dominant at x = and we can therefore also set m? = 0. Thus the operator 
of interest is 

if = -dl + Kflf.Kx) 
VeE/{x) = s{s - l)x-2 + , (3.26) 

It is now straightforward to determine for which values of (p, q) the classical singularities 
at X = become regular or remain singular when probed by scalar waves. First of all, we 
will show that we can reduce the analysis to the case £ = 0: 
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For q < p+2, the first term in the potential is dominant and independent of i. Thus 
A is essentially self-adjoint Ai=q is essentially self-adjoint. As we know from ( 3.1S| ), 
this condition is satisfied iff s{s — 1) > 3/4. 

For q > p + 2, the operators Ai for £ 7^ are essentially self- adjoint by the criterion 
(|3.19|). Thus A is essentially self-adjoint iff ^£=0 is. 



• In the borderline case q = p + 2, for ^ 7^ we have 

^/O 1) > 3/4 (3.27) 

(with equality only for s = 1/2 and i = d = 1). Even in this case, therefore, all the 
Ai with i ^ are essentially self-adjoint and only A£=q needs to be examined. 

We can thus conclude that the operator A is essentially self-adjoint iff s(s — 1) > 3/4 and 

that, in view of ( p. 25 ), it fails to be essentially self-adjoint for 

1 3 2 6 , , 

A not e.s.a. <^ --^ < s < - <^ -- < q < - . (3.28) 

(Jj (Jj 

3.4 The Significance of the (Strict) Dominant Energy Condition 

While this has been rather straightforward, one of the virtues of the present approach, 
based on using a class of metrics appropriate for a generic singularity of power-law type, is 
that it allows us to draw a general conclusion regarding the relation between the Horowitz- 
Marolf criterion and properties of the matter (stress-energy) content of the space-time near 
the singularity. 

Indeed, as we will now show, whenever the matter content of the near-singularity 
space-time is sufficiently "positive" (in the sense of the strict DEC, as it turns out), the 
space-time remains singular according to the Horowitz-Marolf criterion, i.e. when probed 
with scalar waves. 



We can deduce from ( 2.29 ) that metrics with timelike power-law singularities satisfying 
the strict DEC lie in a bounded region inside the strip 0<g<2/(i[||]. Indeed, for g' < p+2 
only the second terms in ( p. 29 ) are relevant, and one finds 

p-Pr = \dq{2 - dq) x~(P+2) ^ (3 29) 

Thus one has ^ 

p-Pr>0 <^ 0<q< - . (3.30) 

d 

In particular, therefore, it follows from ( |3.28D that for such metrics the operator A is not 
essentially self- adjoint and we can draw the general conclusion that 

metrics with timelike singularities of power-law type satisfying the strict Dom- 
inant Energy Condition remain singular when probed with scalar waves. 

Even though metrics with q = 2/d, say, like negative mass Schwarzschild, still satisfy 
the bound ( p. 28 ), thus remain singular while obeying an extremal equation of state, we 
cannot strengthen the above statement to include general metrics with extremal equations 
of state. This can be seen e.g. from examples in [18] and is due to the fact that extremal 
metrics can also be found elsewhere in the {p, g)-plane, in particular in the region q > p + 2, 
while violating the bound ( p. 28 ). 
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3.5 The Friedrichs Extension and "Hospitable" Singularities 

In the previous section we have discussed self-adjoint extensions of (the spatial part A 
of) the Klein-Gordon operator. We have not discussed, however, what these self-adjoint 
extensions imply about the behaviour of the allowed scalar fields (p (those in the domain 
of the self-adjoint extension of A) near the singularity at x = 0. 

It is certainly possible that self-adjointness can be achieved by allowing only scalar 
fields that vanish at the singularity. In some sense, then, the singularity remains excluded 
from the space-time and is not probed directly by the scalar field (p. We will see that this is 
indeed what happens in (in a precise sense) one half of the cases in which there is a unique 
self-adjoint extension. 

However, it is a priori also possible (and perhaps more interesting) to have a well- 
defined time-evolution (which we take to mean "defined by some self-adjoint extension" 
p^) with scalar fields that are permitted to be non-zero at the singularity. In that case, the 
singularity would be probed more directly by the scalar field, and one might then perhaps 
like to define a classical singularity to be "hospitable" (for a scalar field) , if there is a self- 
adjoint extension which allows the scalar fields to take non-zero values at the locus of the 
singularity. We will see that this possibility is indeed realised as well, not only for the other 
half of the essentially self-adjoint cases, but also for e.g. the Friedrichs extension Ap the 
operator ^ in a certain range of parameters for which A is not essentially self-adjoint. 

To address these issues, we need to determine the domain of definition of the relevant 
self-adjoint extension of = —d^ + cx~^ for c = s{s — 1) G [—1/4, oo). For Aq essentially 
self-adjoint, i.e. c > 3/4, this can be done by explicitly determining the domain of the 
closure of the operator Aq . While we have done this (see also ||3^ ) , alternatively, for all 
c > —1/4, one can determine the domain of the Friedrichs extension Ap of Aq, constructed 
from the closure of the associated quadratic form. For c > 3/4, such that is essentially 
self-adjoint, its unique self-adjoint extension of course agrees with the Friedrichs extension. 
Precisely this problem has been addressed and solved in |Q, and instead of reinventing 
the wheel here we can draw on the results of that reference to analyse the issue at hand. 

The main result of [^] of interest to us is their Theorem 6.4. Applied to the operator 
Aq, this theorem^ states that the domain of the Friedrichs extension Ap of is 

D{Af) = {/ G L2(0, oo) : /(O) = 0, / G A{0, oo), d^f G L\0, oo), 

x-^f G L2(0, oo), {-dl + cx'^)f G L2(0, oo)} (3.31) 



where ^(0, oo) denotes the space of absolutely continuous functions. In |33], this result 
was established for c > 0. As far as we can see, this result is correct, as it stands, also for 
—1/4 < c < 0. We will comment on the special case c = —1/4 below. 

We will now extract from this result some restrictions on the behaviour of / near x = 
(assuming that we can model the leading behaviour of / as a; — > by some power of x): 



^Actually, in ^] a more general operator, including in particular a non-zero harmonic oscillator term 
Bx'^, was studied. However, this term serves only to regularise the wave functions at infinity. Since we are 
concerned with the behaviour at a; = 0, this term is of no consequence for the present considerations. 
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1. From the condition x £ we learn that /(x) ~ a;2+'^ for some e > 0. Then the 
conditions /(O) = and dxf G are also satisfied. 

2. The remaining condition (—9^ + cx~'^)f G can be satisfied in one of two ways. 
Either both terms separately are in or / lies in the kernel of the operator (as 
X 0). In the former case, we find the condition f{x) ~ xa"*"^ with e > 0. In the 
latter case, since the two functions in the kernel are x^ and x^~^, with (as usual) 
c = s{s — 1), we now need to distinguish several cases: 



(a) c > 3/4: this means that s > 3/2 or s < —1/2. The solution x^ with s > 3/2, 
i.e. f{x) ~ X2"'"'^, yields nothing new. The solution x^~^ with s > 3/2 (or, 
equivalently, the solution x^ with s < —1/2) is ruled out by condition 1. 

(b) c = 3/4: this means that s = 3/2 or s = —1/2. In this case, we can allow x^^"^ 
and thus relax the domain to include functions /(x) ~ X2+'^, now with e > 0. 



(c) -1/4 < c < 3/4: thus -1/2 < s < 3/2 and s / 1/2. Thus the solution x^ is 
adjoined to the functions {x2+'^} for s > 1/2, and the solution x^~^ for s < 1/2. 



It remains to discuss the special value c = —1/4 or s = 1/2 which is not covered by 
the formulation of the domain in (3.31). This is the minimal allowed value of interest to 



us (c = s{s — 1) with s real), and also the minimal value for which the operator remains 
positive (and thus has a Friedrichs extension). In this case, the two solutions are x* = xz 
and X2 logx, and we checked that, as expected, the domain of the Friedrichs extension 
includes x^^'^. This can also be deduced e.g. from |34], which moreover illustrates nicely 
some of the weirdness of non- Friedrichs extensions. 



The above discussion shows that the two definitions ( p. 14 ) and ( 3.16| ) of the energy, a 



priori differing by boundary terms due to the integration by parts, agree for the Friedrichs 
extension for c > —1/4 and differ only by a finite term for c = —1/4. The issue of boundary 
terms for more general domains is discussed in [^]. 

Returning to the original question of determining the behaviour of the allowed scalar 
fields in the domain of the self-adjoint extension of the spatial part A of the Klein-Gordon 
operator, we need to now undo the transformation (p ^ (j) from the initial scalar fields to 
the half-densities (j) that we performed in section 2.1 to put A into the form of a standard 
Schrodinger operator. 

This transformation back from (j) to (p is accomplished by multiplication by x"''. Now 
the upshot of the above discussion is that the lowest power of x appearing in the domain 
of Ap is 

xi+^ for s > 3/2 or s < -1/2 

x^ for 1/2 < s < 3/2 (3.32) 
x^-" for - 1/2 < s < 1/2. 



~ < 



Evidently these functions are, in particular, positive powers of x. Thus they, and therefore 
all functions in the domain, tend to zero for x ^ 0, consistent with the condition /(O) = 
in ( |3.31| ). However this is not necessarily true for the transformed functions, for which one 



-18- 



has (5 = {5(s) > is a positive real number depending on s) 




x-^ for s > 3/2 

for 1/2 < s < 3/2 

^ for - 1/2 < s < 1/2 

x'^+^ for s < -1/2 



(3.33) 



'mm 



'mm 



1— 2s ™(5 



X 




The final result is the simple statement that a (p in the domain of the Priedrichs extension 
Ap of A necessarily goes to zero for s < 1/2, (p can be non-zero (but remains bounded) for 
1/2 < s < 3/2, and can become increasingly singular for large s > 3/2. 

Note that this statement is not invariant under s — > 1 — s. Indeed, while the operator 
— (9^-|-s(s — l)x~^ has this invariance, and therefore also statements about its essential self- 
adjointness are symmetric under s — > 1 — s (as we have seen), the unitary transformation 
between and 4> depends linearly on s and thus leads to a behaviour of the original scalar 
fields (j) that does not have this symmetry. 

Once again we find a pleasing relation with the DEC, since the watershed happens 
exactly ats = l/2<^q = 2/d which, as we have seen, corresponds to p = Pr- Timelike 
singularities satisfying the strict DEC have < q < 2/d ( p.30| ), thus < s < 1/2. Moreover, 
metrics with s < —1/2 have a unique self-adjoint extension (c > 3/4), thus are regular in 
the Horowitz-Marolf sense, but are not "hospitable" in the sense described above, while 
those with s > 3/2 are. 
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